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Transport through nano-channels plays an important role in many biological processes 
and industrial applications. Gaining insights into the functioning of biological transport 
I/"") ' processes and the design of man-made nano-devices requires an understanding of the basic 

physics of such transport. A simple exclusion process has proven to be very useful in ex- 
plaining the properties of several artificial and biological nano-channels. It is particularly 
qj • useful for modeling the influence of inter-particle interactions on transport characteristics. 

In this paper, we explore several models of the exclusion process using a mean field approach 

and computer simulations. We examine the effects of crowding inside the channel and its 
^(— > . 

t/3 , immediate vicinity on the mean flux and the transport times of single molecules. Finally, we 

+-* ■ 

discuss the robustness of the theory's predictions with respect to the crucial characteristics 

of the hindered diffusion in nano-channels that need to be included in the model. 
Ch ' PACS numbers: 87.10.Ca, 87.10. Mn, 87.85.Rs 

o 
o 

INTRODUCTION 

> 

m 

r**** ■ Molecular transport at the nano-scale lies at the core of many biological processes and novel 

technological applications [l|-|4|. Biological examples of such transport include numerous molec- 

f-*>) ■ ular transporters, such as porins, ion channels, various secretion systems, nucleo-cytoplasmic ex- 

change, and many other cellular channels and transporters [5l4l3]. These biological devices trans- 
port molecules into and out of the cell, as well as between different cellular compartments. This 
transport is highly regulated and selective. Such transporters function as 'molecular filters' that 
transport only their cognate molecules while efficiently filtering out everything else. 

Nano-channels also play an increasingly important role in technological and engineering appli- 
cations. Examples include nano-fluidics, novel biosensors, nano-molecular sorters, nano-pore based 

Qfl.DT 



2l|. The artificial channe 
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s can serve 



22]. 



DNA sequencing and molecular separation processes 

as model systems to mimic and explain the principles of similar biological channels 

Modern experimental techniques provide a wealth of information about the function of such 
devices and an opportunity to directly test the theoretical concepts. Several experimental methods 



are common in the study of the transport through nano-channels. On the macroscopic scale, 



measurements of the static and dynamic be 



with a variety of methods 



,11 Is 



12, 
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ravior of fluxes through the channel are well established 



23H25I] . On a microscopic scale, single molecule tracking 



26M29I]. The various 



techniques can measure directly the transport times and their distributions 
techniques provide information about different aspects of the transport process. 

The mechanisms of the function of such channels are still not fully understood. In order to 
fully realize the potential of nano-channel based technologies and in order to understand the bio- 
logical transport properties, one needs to understand the basic physics of transport through nano- 
channels-from macroscopic fluxes down to single molecules. Despite their apparent biological and 
technological complexity and diversity, it appears that the function of many such transporters can 
be understood from basic principles within simplified models that use a small set of coarse-grained 
parameters and capture the essentials of the problem. These models are commonly based on non- 
equilibrium statistical mechanics |iI.|20U3Gh46], Analytical results obtained from simplified models 
provide important insights into the basic principles of the system and enable one to investigate 
separately the effects of different factors in idealized situations, un-confounded by other factors. 
In addition, such models guide the creation of new artificial channels and can direct one how to 
manipulate the channel properties in order to achieve a desired functionality. When necessary, the 
simple models can be generalized to include more details, such as charge and avidity effects 46M49I]. 
and they may be complemented and compared with molecular dynamics, Brownian dynamics and 



other simulation techniques 
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In general, the choice of a model hinges on the balance between computational feasibility and 
conceptual transparency on one hand, and the ability to capture all the essential richness of the 
experimental setup on the other. One important feature that has to be captured in a model of 
transport through nano-channels is the non-equilibrium crowding and interactions between the 
molecules inside the channel. In particular, the exclusion process model of the confinement and 
coarse-grained interactions of the particles with the channel has proved to be a good starting point 

m 



36, 



52 



55l |. It has been successfully applied to the explanation of various processes, such as 



36 



55] 



re nuclear pore complex 
42H of the hindered diffusion in nano-channels 



the transport of DNA through nano-channels and t 
In this paper, we expand our previous models [36 
as an exclusion process to include the channel asymmetry and non-uniformity, in order to capture 
more realistically the diffusion and binding processes near the channel entrance. The analytical 
results are validated against kinetic Monte Carlo simulations. Finally, we discuss what aspects of 
the model are sensitive to the inclusion of these features and the applicability of the model to the 



analysis of experimental data. 

II. THE MODEL 

In the exclusion process model, the channel is represented as a sequence of N 'sites'. It is 
important to emphasize that these 'sites' do not represent the actual binding sites of the molecules 
inside the channel, but rather serve as a convenient computational tool to model the hindered 
diffusion. The particles enter from the left with constant flux J and can hop to the adjacent 
sites if their occupancy is less than the maximal occupancy m, which models the crowding and 
competition for space. The interactions of the particles with the channel are represented by the 
rates of hopping through and the exit from the channel. This dynamics is illustrated in Fig. [1] 
The model allows the calculation of both the macroscopic properties as the fluxes and transport 
probabilities through the channel as well as properties which characterize the single molecules being 
transported as the average translocation/return times and their distributions. 
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FIG. 1: Schematic illustration of the mapping of transport through a nano-channel onto an exclusion process. 
The rates tl and tr represent the diffusion from the immediate vicinity of the channel entrance towards 
the bulk. r2i and ri2 represent the exit/entrance rates and the dynamics in the immediate vicinity of the 
channel. The rate r represents the diffusion inside the channel. The particles enter at site 1 and hop between 
adjacent sites if their occupancy is less than the maximal occupancy m, which models the steric hindrance. 



III. UNIFORM ASYMMETRIC CHANNEL 

As a simplest approximation, we first discuss a uniform channel, where the dynamics at the 
immediate vicinity of the channel is not considered explicitly (ri2 = T21 = r, see Fig. [1]); we discuss 
the implications of the different choices of the rates in Section IIVI 



1. Single Particle 

We first discuss the transport of a single particle through the channel, in the absence of the 
impinging flux J . The particle starts at the 'entrance' site 1 and can hop to the right to site 2, 
with rate r, or hop to the left, out of the channel, with rate rx,. From any internal site 1 < i < N, 
the particle can hop either left or right with equal probability, with rate r (for each direction). 
From the last site N, the particle can hop either to the right, out of the channel, with rate r R or 
to the left, to site N — 1, with the rate r. The strength of the interaction of the particle with the 
channel is modeled by the choice of ri and r R . It is important to note that particles can hop out 
of the channel in both directions and do not always exit on the right. 

The motion of the particle through the channel is described by the time dependent vector of 
probabilities pi (t) to be at a particular site i: \p(t)) = (pi(t), ...pi(t)...pN(t)), so that (i\p(t)) =Pi(t). 



The Master equation for t 
of the channel ends, is 



30, 



le probability vector \p(t)), describing the hopping through and out 



32 



40. 



m 



±\p(t))=M-\p(t)), (1) 

where the matrix M has the following elements: 

M iti = -2r, M M± i = r for 1 < i < N (2) 

Mi,i = -r - r L , M N>N = -r - r R , M 1>2 = M NjN -\ = r. 

oo 

It is useful to introduce the Laplace transform of the probabilities vector \p(s)} = / e~ st \p(t))dt. 

o 
In Laplace space, the Master equation ([1]) becomes 

s\p( S ))-\p(t = 0)) = M-\p(s)). (3) 

The solution to this equation is \p(s)) = (si — M)~ l ■ \p(t = 0)). Or, in explicit form: 

^^(A'-^r' f + ' + ^-'- A-,), (4) 

where 



and 



2r + s ± Vs 2 + Ars 
Ai,2 = Yr , (5) 



( \i\ -r+ (r + r R + s) Ai 



1/A = - - (r + r L + s) Ai + rXf + ^ , , , . (r + r R + s- r\ 2 ) . (6) 

X 2 J -r+(r + r R + s)\ 2 ) 



The expressions above for the Laplace transform enable one to calculate various quantities of 
interest, such as the transport probabilities and the moments of the transport times. 

The instantaneous probability flux to the right out of the channel is r R pN(t). Thus, the proba- 
bility that the particle exited the channel from the right side by time t is simply the accumulated 
probability flux P^ = L r R p]\r(t')dt'. This is what is often measured in bulk flux or single molecule 



experiments 26]. 



The total probability to exit to the right is 

P^ = P™ = r R (N\e Mt \l)dt' = r R p N (s = 0). (7) 

Jo 

Using the expression for pi(s), we get the translocation probability as 



p rrR fo\ 

~* (N-l) r L r R + r(r L + r R ) ' 



40 



4J. 



As expected, for n, = r R (i.e., for a symmetric channel), this result reduces to known results 
In the limit r R 3> tl and fixed r, P_> —> 1 and when tl 3> r R and r is fixed, P_> = r R /ri — > 0, 
as expected. In the limit of fast diffusion inside the channel, when both rL,r R <C r, one gets 
P_s. = r R j (r R + tl)', that is, the channel essentially behaves as a one-site channel 56J]. Note that 



for an asymmetric channel the translocation probability varies between and 1, unlike the case of 
a symmetric channel where the upper limit for the translocation probability is 1/2. 

We now turn to the calculation of the transport times. The probability density of the exit 
times to the right /-►(£) is the derivative of the cumulative distribution P^., divided by the total 
probability to exit to the right. In Laplace space, this is expressed as /->(s) = - r ^ Pl (s=o) ' Therefore, 
the mean first passage time to exit to the right is 

- 1 dp N (s) 



Pn{s = 0) ds 
Similarly, the mean first passage time to exit to the left is 

t^ 1 dpi(s) 



s=0 



s=0 



(9) 



(10) 



pi(s = 0) ds 
The mean exit time from any of the ends is 

T = -(r L d M^ + r R d M^\ . (11) 

\ ds ds J s=o 

Using the explicit solution for the probabilities, from Eqs. (J4j) ,(j5j) , (JHJ) , we get: 

- N ((N 2 - W + 2) r L r R + 3(JV - l)r(r L + r R ) + 6r 2 ) 

6r((N - l)r L r R + r(r L + r R )) 
_ N ( (2N 2 -3N + l)r 2 R + 6(iV - l)rr R + 6r 2 ) 

<~ ~ 6((JV - \)r R + r)({N - l)r L r R + r(r L + r R )) 
T - N 2 + (N-l)r R /r 

2r(N-l)r L r R /r 2 + (r L + r R )/r' l ' 



Similar results were obtained in [531. 
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FIG. 2: Left panel: Average escape times from a channel versus the interaction energy E = — ln(r^) for 
Tl = 0.01. The solid lines correspond to the average forward escape time and the dashed lines to the average 
backward escape times. The black lines correspond to J = 0, namely the average escape times of a single 
particle and the blue lines correspond to J/r = 1 where many particles are crowded in the channel. The 
channel length is N — 10. Right panel: Comparison with the Kramcrs-Arrhcnius theory. The solid line is 
the uni-directional escape time from a potential well as a function of the interaction energy E = — lii(rn) 
for tl = and J = 0. The dashed line is the exponential asymptote of the Arrhenius-Boltzmann-Kramers 
type. 



The behavior of the mean transport times is illustrated in Fig. [21 The right panel shows 
the escape time from a potential well with a reflective boundary condition on the left (tl = 
0). In the regime of high interaction energy E = — ln(rR), the dependence of the escape time 
asymptotically reaches the exponential exp(— E) in accord with the heuristic expectation of the 



Arrhenius-Boltzmann-Kramers' type 



31 



32|. By contrast, when there is a finite escape rate to 



the left (i"l 7^ 0), the left panel in Fig. [2] shows that the directional average escape time to the 
right does not follow the exponential asymptote as a function of the binding energy E, but rather 
saturates to a constant, determined by the value of the exit rate to the left, ri- 



A. Non-equilibrium flux through the channel: crowding 

1. Steady state flux 

When a finite flux J impinges onto the channel entrance, a non-equilibrium steady state is 
established, so that at any time there can be many particles in the channel that might interfere 
with each other's passage. The particles in the channel obey the same kinetics as described in the 
previous section, with the condition that a site can contain up to a maximal number of particles 
m. One can describe the system in terms of the average site occupancies Uj. After the relaxation 
of the initial condition, a steady-state density profile \n) ss is established. For the case of a uniform 
channel 

n ss = J(r R (N-i)+r) 

J((N - l)r R + r)/m + (N - l)r L r R + r(r L + r R )' 

The average exit flux to the right is </_► = r R n s ^, which yields for the probability of an individual 
particle within this flux to exit to the right 

P_> = J -± = ^ . (14) 

J(l — n\ s /m) (N — 1) ri 1 r R + r (n + r R ) 

Note that the translocation probability of individual particles is the same as in the single-particle 
case (at least for uniform channels), even though the particles are interfering with each other's pas- 
sage through the channel. Another characteristic of the transport is the efficiency which measures 
the fraction of the impinging flux which is being transported through the channel: 

Eff = — = -T, v s-, T ~f- s 1 7- (15) 

J J{{N - l)r R + r)/m + (N - l)r L r R + r(r L + r R ) 



Similar results have been obtained in 



34j, 



42 



54|, 



58]. 



Unlike in a symmetric channel, the exit rates at both ends can vary independently, as they can 
be determined by different physical mechanisms. Moreover, the crowding can also influence the 
diffusion rate r inside the channel (for instance, it has been observed in several systems that the 

n 

on and off rates of molecular binding can change due to crowding [16[ ) . These and other effects 
lead to a rich behavior of the steady state flux through the channel. 

The general qualitative features of the flux behavior are illustrated in Figs. I3|41 In particular, 
they show that in the case when the exit rates vary independently from each other, the transmitted 
flux is a monotonic function of the rates (Fig. [3]) . Only in the case when the values of the exit rates 
are connected (of which a symmetric channel is a special case) via a common physical mechanism, 
the flux (or the efficiency) exhibits optimum with respect to the exit rates (Fig. [3]). Thus, varying 



the relative strength of trapping at the channel ends, it is possible to tune the selectivity properties 
of the channel, which has important implications for the design of artificial nano-channels. 
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FIG. 3: Transmitted flux through a channel of length N = 10 versus the escape rates r^, tr. The impinging 
flux is J/r = 0.1 for all lines except the black dashed line in the right panel for which J = 0.1 in the arbitrary 
units used for the other rates (in this line the hopping rate inside the channel r depends on the hopping 
rate outside the channel tr). 



The dependence of the transmitted flux on the impinging flux J is illustrated in Fig. [U The 
left panel of Fig. [5] shows the dependence of the transmitted flux on the impinging flux when the 
diffusion rate inside the channel depends on the impinging flux (the different lines correspond to 
different dependence of the internal diffusion on the impinging flux) and the rates of hopping from 
the channel outside are fixed. When the diffusion rate inside the channel is independent of the 
impinging flux or is a monotonically increasing function of it, the transmitted flux exhibits the usual 
saturation curve. Such a situation can arise when the effective binding energies decrease due to 
inter-particle crowding [16|, |59(. On the other hand, when the diffusion rate inside the channel is a 
monotonically decreasing function of the impinging flux there is an optimal value for the impinging 
flux which maximizes the transmitted flux. Such a situation can arise experimentally when the 
transported molecules change the conformation and structure of the channel constituents [26 1. 

The crowding can also have an effect on the exit rates. In the right panel of Fig. [H we show the 
dependence of the transmitted flux on the impinging flux when the exit rates out of the channel 
directly depend on the impinging flux. Similarly to the case where the flux influences the diffusion 
inside, the transmitted flux exhibits an optimum as a function of the impinging flux both for the 



case in which the exit rates are a monotonically decreasing function of the impinging flux and 
for the case in which they depend quadratically on the impinging flux. When the exit rates are 
independent of or increase with the impinging flux the optimum disappears, and the transmitted 
flux monotonically increases and saturates as a function of the impinging flux. 
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FIG. 4: Transmitted flux through a channel of length N = 10 versus the impinging flux. The kinetic rates 
can depend on the impinging flux; see text. Left panel: tr = 0.1, tl = 0.01. Right panel: r = 1. 



2. Single particle on the background of the steady state flux 



In order to interpret single-molecule tracking experiments in the experimentally relevant 
regimes, we investigate how the transport times of the individual particles are influenced by the 
crowding, when many particles are present in the channel. To the best of our knowledge, no closed 
analytical solution exists in this case. We use the mean field approximation, backed by kinetic 
Monte Carlo simulations. 

To this end, transport of an individual particle can be viewed as occurring on the background 
of the steady state density profile \n) ss . In the mean field (or the effective medium) approximation, 
the probability pi (t) of a particle to be present at a given site is described by the following equations 
4J: 



— pi = rp i -i(l-nt s /m) + rpi + i(l-ni s /m)-rpi(l-nli 1 /m)-rpi(l-nt+ 1 /m), (16) 



10 



with the boundary conditions 
d 



dt 
d_ 

It 



Pi = rp%(l - n\ s /m) - rpi(l - n s 2 s /m) - r L pi, 



(17) 



PN = rp N -x{l - rix/m) - rp N (l - n s ^_Jm) - r R p N . 



Similar to the case of a single particle, this set of equations can be written in a matrix form: 

d 



dt 



\p(t)) = M ss -\p(t)), 



(18) 



with the matrix elements of M ss following from Eqs. (|16ll7p . In Laplace space, this becomes 

\p( S )) = (si-M ss )- 1 -\p(t = 0)). 



(19) 



3. Mean transport times 

Equations (|16p .(|17 p and (|18p can be solved in a manner similar to the single particle case above. 
In particular, the translocation probability turns out to be the same as for a single particle in the 
absence of flux, as has already been established in the previous sections (based on the ratio between 
the transmitted and entering fluxes). After some algebra, one gets for the mean time to exit to the 
left: 



N 



^+4^ + 3(iV-l) 



JP- 



+ ^ R + N 



n^+^-^T^) 



2 J IN 



^)/- 



.(20) 



where i/)(x) = ^r(x); T(x) is a 7-function. A similar expression can be obtained for Tff (not 
shown). The mean time to exit to the right can be calculated from the relation 



p^ri + p^t: 



(21) 



Interestingly, the mean escape also turns out to be unaffected by the crowding, similar to the 



uniform channel case 



42| (see Fig. E]). 



The dependence of the mean escape times on the impinging flux is shown in Fig. [5l Similarly 
to the symmetric channel case, the crowding increases the forward exit time and decreases the 
backward exit time. The mean field approximation is compared with the kinetic Monte Carlo 
simulations (dotted lines in Fig. EJ), which shows that it is an excellent approximation for wide 
channels (m = 3) and is a reasonable approximation also in the case of single file transport (m = 1) 
even for quite high values of the impinging flux J, where the crowding effects are significant. 



11 



Fig. [2] (right panel) illustrates the dependence of the average escape times to the left and to 
the right on the interaction energy. Interestingly, the crowding does not change the qualitative 
behavior of the escape times. 

It is also interesting to note that 



J(l - nl s )T 



i occup ) 



(22) 



where N, 



yi=N 



occup 



Yll=i n i i s the average number of particles in the channel. Heuristically, despite 



the crowding effects, one can think of the channel as a 'box' with 



i occup 



particles that hop out 



with an average rate 1/T so that the total flux out of the channel is «/_» + J<. 



N, 



occup 



/T ss . This 



heuristics is very different from the one applied to the potential well with a uni-directional escape 
(see Kramers' theory). For instance, it is not true that </_>. = N occup /T^.. It is important to note 
that the above mentioned heuristics is only valid when thinking about the average flux; higher 
moments will deviate from that picture since the distribution of exit times is not exponential. 
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FIG. 5: Dependence of the mean exit times on the impinging flux J, normalized to the case with J = 0. 
Blue middle line: the average dwelling time T^/T^. Lower line: the return time T^_/T^. Upper line: 
the forward exit time T^/T^. The corresponding simulation results are plotted with dotted lines. The 
parameters are: N = 6, r = 1, m = 3, r = 1, rx = 0.01, rp = 0.1. 



4- Probability distributions of the transport times 



The single molecule tracking experiments can measure not only the mean transport times, 



but also the transport time distributions 



26. 



28, 



29l | . which can be calculated within the present 



12 



formalism. As mentioned above, the distribution of the exit times to the right /-►(£) is simply 
r RPN(t)/P^, where p?j is the probability that the tagged particle is in site N at time t and the 



30 



31 



, 142], |57|]. Thus, from 



distribution of the exit times to the left is /<-(£) = rLPi{t)/P+ 

equations ([TBI [T71 [T8j) , the probability distribution of the exit times to the right of a particle that 

enters the channel at site 1 is 

U(t) = ^exp(M ss t) . (23) 

JT^jp \ / ivl 

In Laplace space, we can write it as f->(s) = (rji/P^)(N\(sI — M ss ) _1 |p(t = 0)). In the single 
particle limit, J — > 0, when the n ss — > and M ss — > M, it reduces to the expressions of the previous 
section. The probability distributions of the transport times are illustrated in Figs. El [8] and [71 

Both the forward and the backward exit times have an exponential decay tail (as demonstrated 
by the straight lines in Fig. (6]). Note that the decay constant is the same for both forward and 
backward times (as demonstrated by the fact that the lines for the forward and backward exit 
times probability density are parallel in Fig. [6]), because the tails of the distributions are due to 
particles that have been bouncing back and forth inside the channel for a long time and thus have 
lost the memory of the initial condition. 

The distribution of the return time is strongly non-mono-exponential. The decay at very short 
times is much faster than at long times. This can be understood by noting that the decay at 
very short times is due to those returning particles that jump out immediately after entering the 
channel. Therefore, the short time decay rate of the probability density of the left exit times is 
r + 77. It is important to note that even in the limit of simple diffusion, when 77, = tr = r, 



the distribution of the backward exit times is still non-mono-exponential: t 



lis bears important 



28]. 



implications for the interpretation of single molecule tracking experiments [261 . 

Comparison of the mean field probability density with the simulations in Figs. [7J and [8] shows 
that for m = 3 the mean field and simulation curves are almost identical for all values of J. For 
the single file case, m = 1, and large values of J there is a small difference between the mean field 
and the simulation curves due to the fact that the mean field fails to capture correlations between 



successive hops for strictly single file transport in the crowded regime 



42] 



IV. NON-UNIFORM CHANNELS: PARTICLES RETURNING TO THE CHANNEL 

The uniform channel is one of the simplest models that captures the essentials of hindered 
diffusion with interactions through a narrow channel. It has been shown to account well for 
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FIG. 6: The distributions of both the forward and the backward exit times decay exponentially with the 
same constant at long times. The parameters are shown in each panel. In all panels TV = 6, r = 1, and 

Th = Tr = 0.1. 



the effects of the inter-particle interactions on the flux 



361 ]. and some of its more complicated 



predictions in the multi-species case have been verified experimentally 22|, [55(. An interesting and 



useful feature of the uniform channel model is that man y c 
the impinging flux and are un-affected by the crowding [36 



uantities of interest are independent of 



42, 



m- 



However, when thinking of a uniform channel as a model of diffusion through a channel, it is 
apparent that it neglects an important effect: the possibility that a particle can return into the 
channel after it has exited. Moreover, the diffusion rate in the immediate vicinity of the channel 
entrance can differ from both the diffusion rate in the bulk and inside the channel. We now model 
the dynamics near the channel entrance in more detail, by accounting for the possibility that the 
particles can return into the channel after they exited it. This effect is modeled by choosing the 
sites 1 and N to lie 'outside' the channel. This is illustrated in Fig. [TJ Now the rate T2\ represents 
the escape rate from the channel to the medium in the immediate vicinity of the channel end, 
while the rate r\2 represents the return rate into the channel. Accordingly, r^ and r^ represent 
free diffusion from the immediate vicinity of the channel towards the bulk (to the left and right 
correspondingly). As before, r is the diffusion rate inside the channel and J is the flux impinging 
from the bulk. The sites outside the channel (1, AT) are assumed to be free of any restriction on 
the number of particles they can contain. 

The local exit and entry rates r<i\ and ryi contain the potentially complicated dynamics of 
binding/unbinding and local diffusion near the entrance. From the detailed balance condition, 
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FIG. 7: Backward exit time distributions for different values of J. Solid blue line: mean field results. 
Dotted red line: simulations. Upper panels: m=l;Zower panels: m=3. In all panels N = 6, r = 1, and 
r L =r R = 0.1. 



their ratio is equal to the Boltzmann factor of the energy difference E between the inside and the 
outside of the channel: r 2 \jr\ 2 = exp(-E). However, their actual values are determined by the 
local diffusion coefficient, which in turn depends on the details of the energy profile of the molecular 
interactions near the channel entrance. 

The mean field equations for the steady state densities in this case are [32j, 
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52, 



53] 



= -2rnf + r (nf^ + nfi x ) for 3 < i < N - 2, 

and for the boundary sites: 

J = r L n{ s + r i2 nf (1 - nf/m) - r 2 in s 2 s 

= — r 2 in| s - rn| s (l - nf/m) + ri 2 nf*(l — n s 2 s fm) + rnf(l — nf/m) 

= -r 2i n^_ 1 - rn5J_!(l - nj$_ 2 /m) + n 2 n$(l - nf-i/m) + rn^_ 2 (l - nf^/m) 



(24) 







-ran" - r 13 n#(l - n'$_ x /m) + r 2 m s ^_ l . 



(25) 



These can be easily solved by noticing that the mean field equations for the densities at the 
internal sites are identical to the equations describing the steady state population in uniform 
channels: nf s = a + bi for 2 < i < N — 1. The set of equations (|24p can then be solved for the 
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FIG. 8: Forward exit time distributions for different values of J. Upper panels: m=l; Lower panels: m=3. 
In all panels TV = 6, r = I, and r^ = t\r = 0.1. 



variables a,b,nf s and njy. The explicit solution has been obtained in Mathematica but since it is 
cumbersome it will be omitted here. 

The flux to the right through the channel is given by J_> = rRn^| and the transport efficiency 
is given by </_>./</, which can be readily calculated from Eqs. (|24j) . Not all the particles in the 
impinging flux are able to enter the channel, because the channel entrance (site 2) can be occupied. 
Thus, the transport probability of a particle to traverse the channel to the right once it has entered 
(into the site 2) is J^./J- m , where J m = J I 2 f^^i/L is the flux that actually enters the channel. 

Figures [9] and [10] show the transport efficiency and the translocation probability as a function 
of the rate t<ix for fixed ratio rixl^xi- That is, in those figures the local diffusion coefficient at the 
channel entrance/exit varies, while binding energy stays constant. For wide channels (maximal 
occupancy m = 3), there is an excellent agreement between the mean field calculations and the 
results of simulations of the exact dynamics. For m = 1 (single file dynamics) the mean field 
still describes the efficiency reasonably well. The transport probability decreases with r2i, because 
larger t<i\ implies an increased chance of hopping out of the channel. Somewhat counter-intuitively, 
the transport efficiency increases with ri\. This is due to the fact that for fixed ratio r%\jr\ii the 
larger r%x implies larger ri2, which acts to keep the particle near the channel entrance instead 
of diffusing away to the bulk. These results indicate a possibility of an additional 'knob' for the 
control of transport properties of nano-channels. 

Figures ITU and [T2l show the efficiency and the translocation probability versus TiXt but in the 
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FIG. 9: The transport efficiency J^/J versus the local diffusion coefficient at the channel ends for a fixed 
ratio T2i/t\2 = 0.1. The blue lines show the analytical mean field results while the red lines show the 
simulations of the exact dynamics. The left panel shows the case where the maximal occupancy of each 
site is to = 3 and the right panel shows the case of single file channel to = 1. In all panels J = 0.1, 
r = r^ = rn = 1, and TV = 10. 
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FIG. 10: The translocation probability P^ versus the local diffusion coefficient at the channel ends for a 
fixed ratio ri\jr\i = 0.1. The blue lines show the analytical mean field results while the red lines show 
the simulations of the exact dynamics. The left panel shows the case where the maximal occupancy of 
each site is to = 3 and the right panel shows the case of single file channel m = 1. In all panels J = 0.1, 
r = Ti, = r_n = 1, and N = 10. 
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case where r\i is fixed, that is as a function of the binding energy at the channel entrance. As 
above, the mean field approximation provides a good fit to the results of the simulations. However, 
unlike in Fig. [9l the transport efficiency is a non-monotonic function of ri\ and has an optimum 
at a certain value of r2i ■ In thi s respect, the non-uniform channel behaves as a uniform case that 
has been studied before 
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FIG. 11: The transport efficiency J->/J versus the binding energy at the channel ends. For these plots we 
have used the following parameters: J = 0.1, r = r^ = tr ■ = r\i = 1 and N = 10. The left panel shows the 
case where the maximal occupancy of each site is to = 3 and the right panel shows the case of single file 
channel m = 1. The blue lines show the analytical mean field results while the red lines show the simulations 
of the exact dynamics. 



Diffusion outside: choice of the model 



To what extent can the effect of the returning particles be neglected when analyzing the exper- 
imental data? Clearly, the non-uniform model, in which the sites outside the channel are explicitly 
considered, is more complicated and involves more parameters than the uniform channel model. It 
is important to know whether bulk measurements of the flux can differentiate between the two mod- 
els, necessitating the explicit inclusion of the particle return for interpretation of the experimental 
data. 

Here we start investigating this question by comparing the mean field exit flux through a 
symmetric non-uniform channel with length N and rates r\% , r-i\ , r r = r^^r with the exit flux 
through a uniform channel with the length N and rates r , r. Qualitatively, both fluxes exhibit a 
saturating Michaelis-Menten type behavior as a function of the impinging flux J. It turns out that 
with the appropriate choice of the effective parameters N and r the behavior of the fluxes is also 
quantitatively similar. 
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FIG. 12: The translocation probability versus the binding energy at the channel ends. The parameters are 
the same as in Fig. 111! The left panel shows the case where the maximal occupancy of each site is m = 3 
and the right panel shows the case of single file channel m = 1. The blue lines show the analytical mean 
field results while the red lines show the simulations of the exact dynamics. 
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FIG. 13: The output flux through a non-uniform (solid blue line) and an approximating uniform channel 
(dashed red line) versus the impinging flux J. The parameters are r = ?\r = Tl = 1, ^12 = 2, r%i = 0.5 and 
2V = 10. The fitted parameters are N = 7.8 and r = 0.3, close to the naively expected values (see text for 
details). 
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This result is illustrated in Fig. [13j which shows the exit fluxes versus the impinging flux for the 
nonuniform model and the corresponding uniform model. The effective length N and the effective 
exit rate r for the uniform channel were obtained by requiring that the exit flux at the saturation 
limit J — > oo and the slope of the exit flux with respect to J at J = are equal for both models; 
these are two parameters that completely determine the shape of the Michaelis-Menten type curves. 
In Fig. [T3l the parameters for the non-uniform model are: r = tr = rj_, = 1, r\2 = 2, r<i\ = 0.5 
and N = 10. The fit of the uniform model yields the following parameters: N = 7.8 and r = 0.3, 
which are not too far from the effective parameters one would obtain by naively removing the 
sites outside the channel (in the equivalent uniform model) and setting the ratio r /r = T2i/r\2- 
The fitted curve is very close to the curve obtained from the non-uniform model (the deviation is 
much smaller than the experimental precision in the experiments that we are familiar with). The 
agreement between the two models is rather general and does not depend on a fortuitous choice 
of parameters (data not shown), although its quality decreases when ri2,i"2i *C r. Therefore, in 
general, in order to distinguish between the different models, additional experiments at the single 
molecule level are needed. 

It is important to note that the two lines coincide at the limits J — > and J — > oo (these are 
the fitting requirements) but still deviate in the intermediate regime. This deviation is a clear 
indication of a non Michaelis-Menten behavior of the output current for non-uniform channels. 
These effects lie outside of the scope of the present work and will be studied elsewhere. Finally, it 
is important to compare the models in the experimentally relevant regime. We have recently shown 
that the exclusion process model describes well the observed behavior of the flux of single-stranded 
DNA through nano-channels functionalized with DNA hairpins [la, [3(| . For the parameter values 
arising from those experiments, the two models are essentially indistinguishable (data not shown) 
for the same values of the physical parameters. 

V. SUMMARY 

In this paper we have investigated the behavior of the exclusion process models of non- 
equilibrium transport through nano-channels, using a mean field (effective medium) approach and 
kinetic Monte Carlo simulations. We have investigated the dependence of the static and dynamic 
characteristics of the transport on various parameters such as the binding energy in the channel, 
the diffusion rate inside and outside the channel and the flux through the channel in a uniform 
asymmetric channel. In particular, we have shown that the analytical mean field approximation 
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captures well the properties of the transport both at the macroscopic level (fluxes, probabilities) 
and the single molecule level (transport times and their distributions). The results provide addi- 
tional insights into how the properties of the artificial nano-channels can be manipulated to achieve 
a desired functionality. 

We have also investigated an extended model which explicitly takes into account the dynamics 
of particles' exit and return into the channel in the immediate vicinity of its entrance. For a wide 
range of parameters, in terms of describing the behavior of the steady state properties, such as 
fluxes and transport probabilities, inclusion of this effect does not result in a difference that could 
be distinguishable on the basis of the measurement of the flux alone. However, the effects of local 
non-uniformity near the channel entrance might be observable in single-molecule experiments. 

To summarize, simple physical models of the diffusive transport through nano-channels pro- 
vide important insights into the mechanisms of their function. Even the simple exclusion models 
described in this paper have a rich behavior that appears to be capable of grasping the essential 
physics in several experimental systems. Nevertheless, the question of choice of the pertinent pa- 
rameters that have to be included in the model for meaningful interpretation of the experiments 
still remains open and requires future works, both theoretical and experimental, such as single 
molecule tracking measurements. 
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